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ABSTRACT 

This study was conducted to investigate the effects of several parameters in the Navier - Stokes equations 
describing the nanofluid flow subjected to different electric fields. The boundary layer of SWCNTs - water flow over a 
stretched surface ofdifferent thicknesses in the presence ofweak and strong electric fields was analyzed. The governing 
equations were essentially a set of ordinary differential equations (ODEs), solved using the R.K Fehlberg method with 
shooting technique (via MAPLE 18 software). The positive correlation between electric field and temperature was rather 
significant. The mass transfer rate increased with respect to most of the working parameters, indicating that the mass 
transfer is dominated in the presence of electric field with stretching sheet in different variables. 

KEYWORDS: MHD Carreaunanofluid Flow, Single Walled Carbon Nanotube, The Thermal Conductivity, Electric Field, 
Stretching Sheet With Deferent Variables & Shooting with Rungekutta Method 

Receiyed: Dec 23, 2018; Accepted: Jan 13, 2019; Published: Feb 27, 2019; Paper Id.: IJMPERDAPR201934 

NOMENCLATURE 


1 -m 

B(x) -B 0 (x + b) . M a g ne ti c f ac tor 
° : Electrical conductivity 

1 -m 

E(x) -E 0 (x + b) . e i ec t r i ca i fi e ld factor 
u,v : x- and y- velocity 

Q 

w : Wall concentration 
q 

00 : Ambient concentration of the nanoUuid 
T : Nanolluid temperature 
C :Nanofluid concentration 

^ : Nanolluid effective density 
cc f 

n ' : NanoAuid thermal diffusivity 


www.tjprc.org 


SCOPUS Indexed Journal 


editor@tjprc. org 


Original Article 








350 


Jaafar Abdul Abbas Abbood Al-Nasrawi, Rozaini Bin Roslan, 
R. Kandasamy, Inas R. Ali & Radiah Mohamad 


^ n f : NanoAuid dynamic yiscosity 

[pc p ) 

nf : NanoAuid heat capacitance 

1. INTRODUCTION 


The thermal conductivity of nano refrigerant increases significantly with respect to its volume fraction, [1]. In this 
research, the volume space formation of nanoparticle cluster in fluid was fabricated and the impact of the thickness of the 
intake sheet of nanoparticles was studied [2]. The formation of the single walled carbon nanotube could be regarded as a 
single smoothie, curved up graphene sheet. The graphene sheet contains sextuple grids of carbon atoms. Single walled 
carbon nanotube (CNT) contains of a single graphene sheet with exemplary diameter 1.0 -1.5 nm. On the other hand, 
multi-walled (CNT) (MWCNT) contains many graphene sheets with thickness curved up one inside. The typical diameter 
of (CNT) is very much dependent on the size of the impetus particles used to produce the (CNT)[3-7].SWCNTs-water has 
been used in nano biological gadgets and nano mechanical frameworks, [8-12]. They are generally adopted to improve the 
solubility. Meanwhile they can be stabilized by a surfactant within the base fluid to prevent agglomeration and 
sedimentation[13-14]. Several researchers including [15-17] have used nanoAuids to improve heat transfer in thermal 
system of different geometries subjected to different thermal boundary. Some related studies in micro channels have been 
conducted as well[ 18-27]. Thermal radiation energy due to SWCNTs on MHD nanolluid flow in the presence of 
seawater/water: Lie group transformation, [28]. 

The aim of the current work was to study the mathematical scheme used to analyze the Carreau SWCNTs - water 
nanolluid flow due to nonlinear stretching sheet with extensible thickness. However, there have been no controlled studies, 
comparing the differences between the work in (Khan, Mair et. al 2017) and our previous work in simulating Carreau 
SWCNTs - water nanolluid flow due to nonlinear stretching sheet with extensible thickness [29]. 

2. BASIC EQUATIONS 

The aim of this study was to clarify several aspects of the mathematical equation used to study the Carreau 
SWCNTs - water nanoAuid flow due to the nonlinear stretching sheet with extensible thickness. The sheet was developed 

1-771 , 

with velocity U w = U 0 (x + b ) m , where U 0 is the source velocity. B(x + b)~^~ , is the thickness of sheet and ® is the 

dimensionless constant, m is the power law index. It was known earlier that the model must be pleased only for m ^ 1, as 

the model reduces to using flat sheet as for m = 1 . The flexible magnetic field strength ^ was imposed vertically to the 
plate. The electric and magnetic fields were acting perpendicular to the surface. 
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Figure 1: Physical Configuration of SWCNTs - Water NanoAuid System 

Equation (1) is the continuity equation where the current function is used to describe the x- and y- yelocities. 


du _ dv 
dx dy 


du du Mnf 
dx oy p nf 




d 2 u 3(n-l) 2 

df + 2 


du 

dy 


d^u 


\ 


^(e(x)B(x)-B\x)u) 

Pnf 


dT dT 
U dx +V dy 


k nf 

(d 2 T\ 

1 dc ’r | 1 x 

{ s jacarl (A-l/l 

'dT^ 

11 

(P Cp ) nf 

W) 

(PC p )nf (P c p)„f 

^ B ’ nf \dy dyj ] 


11 


+ Ofr-rA 

(. P c p)nf 


— + v— 
dx dy 


= {D, 


B’nf 






{Pr)nf\d 2 T 

T„ |ay 2 


+ K(C„-C) 


The correlated boundary equations for stretching sheet are 

u [x,B(x + by 2 = \J 0 (x + b) m ; v[x + B(x + bf 2 = 0 ; T(x + B(x + b+ 2 ’ ) = T w 

C (x + B(x + = C w ; u(x, 00 ) -> 0 ; T -* T a C -* C^, as y -> 00 


( 2 ) 

(3) 

(4) 


(5) 


Pnf (^ C)Pf + (sP s ’ Pnf 


Pf 

d-O 2 - 5 


(pP)nf = (1 - O(PP)f + C(PP)s . (P C p ) nf = (1 - C)(P c p)f+ C(P c p)s 


01 nf = 


k _4 

(P c p)nf 



\ (k s +(l-l)k f )-(l-mkf-k s ) } 

\ (k s +(i-\)k f )+C(kf-k s ) \ ,<Tnf 


= Cg s + (1 - C)Of (D B ) nf = (1 - C)(d b ) f 


( 6 ) 


The Maxwell design method was used to assign the dynamic thermal or electrical conductivities of liquid-solid 
suspensions. Here, o f and o s are the electric conductivities, k f and k s are the thermal conductivities, C, is the 

nanoparticle volume fraction, p f and p s are the densities, p f is the dynamic viscosity and k nf is the effective thermal 
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conductiyity. Note the Rosseland's dissipation for the radiatiye heat flux is = 4 ^r 4 [3Q]). where&* is the Rosseland 


3 k dy 


mean absorption coefficient and cr is the Stephen-Boltzmann constant. Finally, the termT 4 can be experded from too via 
Taylor's series r 4 = 4rjr - 3ri . 

via similarity conversions: 


¥~~ 


2v f U 0 (x + b)‘ 
m +1 


m -1 


-/W^= 






T -T 


CL-CL 


(7) 


w =-,v =- 

The stream function y/ can be constructed as ^ which adequately attends the Equation (1). 

Depending on equations (6, 7), equation (1) is fairly inlluenced while the equations (2) and (4) incomes the representative 
of three works that will be examined, all of which 


f 

1+ 


\ 

+ G.B 

2m F' 2 FF" 

V 

v 2 

/ 

K m +1 j 


+ G.A.M(E 1 -F' 2 ) = 0 


0” + Vr f .E (f 0’ + (A\)Nb 6'<!>' + (A\)Nt 6' 2 + Ad)= 0 


f + (A\). 


1 Nt 


0" + LeF 6' - Le y(f> 


= 0 


(Al) Nb 
The boundary conditions are 

(1 + m)F(a) = (1 -m) a , F\a) = 1, 0(a) = 1, (j)(a) = 1; F'(oo) = 0, 0( oo) = 0 ,</>(<*>) = 0 


K nf = (k s + (e -1 )(k f )-(e- 1) C(k f - (k s ) K f = (k s + (e - \)(k f ) - £(k f - (k s ) , B = 


p f I K, 


K. 


( 8 ) 

(9) 

( 10 ) 

( 11 ) 

( 12 ) 


(13) 


i-<r+<r 


Ps+P), 

Pf(cp)f 


A = 


1+ 


f \ 

5--1 


f \ 

^ + 2 

a f 

V v j y 


^-i 

. 


,G = {\-Cf- 5 ,A\ = 


i-<r 


Here, a =B 


t/n 


is the related wall thickness parameter and n = K +o( OT+ 0 is the dependent plate surface. 


2v, 


2Vj-(m+l) 

Up on applying transformations, f(a) = f(r/ -a) = f(ij), 0(a) = 0(r/ -a) = 0(r/), and </>(a) = (/)(?] -a) = (/)(r/), 
Equations (8-10) become 


r 


i+ 

V 


^nl w \ r y 


W f 

+ G.B 


/J 


\m +1 j 


+ G.A.M(E 1 -f' 2 )=0 


6" + Pr f .E.(f6' + (A\)Nb 6’<t>' + ( A\)Nt 6' 2 + X6 )= 0 


(14) 

(15) 
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(j)" + (A\). 


1 Nt 

- 0" + LefO' - Le yt/> 

(Al) Nh j 


= 0 


(16) 


1 — 171 

f(a) = a --, f\a) = 1, 6(a) = 1, <t>(a) = 1 ;/'(«=) = 0,0(oo) = 0, *(oo) = 0 

1 + m 


(17) 


„ 2c7 f B 0 2 (jc + fo) 1_m 

Here, M = f 


U+c + b) 3 ™-'- 


, , is the magnetic parameter, t+ = (m + 1)- 

U 0 {m + l)p f 2v 

(jUc 


r 2 is the Weissenberg number, 


Pr P f J K 

p r , -_ is the effective Prandtl number where Pr =- is the Prandtl number, y = - — -is the 

f (*- • ' kf ' U 0 {m + \\x + b) m ~ l 


^+ 4 R 

K k f 3 V 


chemical reaction, Re = 


xU 0 {x) 


Reynolds number, £, 


(pc) (D B ) f (C -C^) . . . ^ 

N h =- --is the Browman motion, N f 

(pc) f Vf 


B 0 u 0 { x +by 

(P C ) S (D T ) f (T w —T„) 


is the electric parameter, 


(pc) f v ^ 


v f 

is the thermophoresis, Le = ^ j is the 


n_ 

Lewis number and K ~ is the solar radiation-conduction. Besides, the physical resources of engineering hits of 


k + 

materials boosting application are skin friction coefficient C f , Nusselt number Nu x and Sherwood number Sh x expressed 
as 


R &\C f = t— 

/ (, _ ^)2.5 


171 + 1 


f{0), Re~ 2 Nu x =- 


k -fL+ 4 -R 


\ k f 


171 + \ 


6'(0 ), Sh x Re , 2 =- 


(1-<T) 


m + 1 


/( 0 ) 


(17) 


3. DISCUSSIONS 


Figure 2 indicates that the magnetic field (M= 0.1,0.5, 0.9) effect with high electric field, the temperature profile 
decreased and such that, the temperature increased with respect to the strength of magnetic field when the electric field is 
weak. Table 2 shows that the rate of skin friction increases when the strengths of both electric and magnetic field are high. 
the heat transfer (—0 (0) ) decreases in weak electric field and strong magnetic field, it is decreased but it is increasing with 
high electric field that means the magnetic field is dominant on the heat transfer. Table 2 shows that the concentration 
decreases with respect to the strength of electric field in the presence of magnetic field. 

Figure 3 shows that when the wall thickness parameter increases, the temperature decreased with high and low 
electric field, and the thermal thickness is large. Table 3 indicates that skin friction increases with respect to the strength of 
the electric field. As shown in Figure 4, in weak electric field, the temperature is independent on the power-law index. 
Nevertheless, the change in temperature is more apparent when the electric field is strong (i.e. n=1.2, 2, 3) power law index 
increased no change on temperature profile with low electric field and there is slightly change when the electric field is 
high such that n=1.2, 2, 3. Table 4 shows that skin friction decreases with respect to the power law index when the electric 
field is weak and when the temperature is low. Meanwhile, the temperature increases when the thermal radiation increases 
low electric field in the same trend with high electric field the temperature increases and the difference between low and 
high electric are slightly decreased. It is interesting to note that the thermal radiation thickness is large indicating that the 
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nanoAuid is dominated by the thermal radiation. 

From Table 5, seemingly the skin friction f (0) is independent on the strength of the electric field. For the heat 
transfer the table shows a decrease with different electric field and an increase for the rate of mass transfer. Figure 6 
indicates that the temperature change is only apparent at high Weissenberg number (i.e. We = 0.5, 1.5, 3.5). Table 4 shows 
that the skin friction f (0) decreases at high Weissenberg number when the electric is low with decrease in 
temperature—0 (0), but increased the rate of mass transfer with different electric field. 

As shown in Figure 7, the temperature drop is more apparent when electric field is strong. However, the 
temperature increases when the Brownian motion is significant (i.e., Nb= 0.01, 0.1, 0.3). Table 7 shows that the drop in the 
rate of skin friction f (0) is more apparent as the electric field is weak. In addition, the heat transfer —0 (O)increasedwith 
low electric and take two behaviour with high electric field. As shown in Table 7, the increase in mass transfer is more 
apparent when the electric field is weak. From Figure 8, the temperature increases with respect to the thermophoresis 
parameter in the deference electric field, is such that the deference between E 1=0.0 and E 1=5.0 are slightly decreased. 
The change in skin friction is significant when the electric field is strong. As reported in table 8, the heat transfer rate is 
rather low when the electric field is weak. Figure 9 shows that the temperature is independent on the Lewis parameter. 
On other hand the deference between the effect on Lewis parameter with low and high electric field is slightly decreased. 

Table 9 indicates that the heat transfer rate—0 (0)is low in strong electric field. The mass transfer rate, however, 
increases with low and high electric field. Figure 10 shows no effect of low and high electric field with increased chemical 
reaction, and there is slight decrease with high electric field. The change in skin frictionf (0), is insignificant is strong 
electric field. Heat transfer drops when E 1=5.0. Finally, mass transfer increases with deferent electric field. Figure 11 
indicates that temperature increases as the heat source increases. On the other hand, the temperature decrease with high 
electric field compare with low electric field. Table 11 shows that heat transfer —0 (0) decreases with deferent electric 
field. Figure 12 shows that, the prandtl parameter trend conversely. It is apparent that Prandtl parameter increases then the 
temperature decreases, independent on the strength of the electric field. Table 12 shows that heat transfer rate would 
increase if the electric field is strong. 


Table 1: Thermo-Physical Properties of the Fluid and Solid Particles 



p(kg/m 3 ) 

C P (J /kgK) 

k(W/mK) 

<7(n _1 m _1 ) 

Pure water 

997.1 

4179 

0.613 

5.5 

SWCNTs-water 

2600 

42.5 

6600 

1.26 


particle shapes 


Sphere, (m = 3.0) 
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SWCNTs-water 

Solidline-E 1 =0.0;Dashline-E 1=5.0 



n 


- M=0.1 - M=0.5 - M=0.9 - M=0.1 

- ~ M=0.5 - M=0.9 _ 


SWCNT s-water 

Solidline-E 1 =0.0;Dashline-E 1=5.0 



n 


- oc=0.1 - oc=0.5 - or=0.9 - o=0.1 

■ ~ o=0.5 - o=0.9 _ 


SWCNT s-water 

Solidline-El =0.0;Dashline-El=5.0 





- n=1.2 -n=2 - n=3 - n=1.2 - n=2 

- n=3 _ 


Figure 2: Magnetic Field (M) 


Figure 3: Wall Thickness (a) Figure 4: Effect of Power Low Index (n) 

on Temperature in Different Electric Fields 


SWCNTs-water 

Solidline-E 1 =0 0;Dashline-E 1=5.0 


SWCNTs-water 

Solidline-El=0.0 ;Dashline-El=5.0 


SWCNTs-water 

Solidline-El=0.0;Dashline-El=5.0 







— — w P= n > wr— i ^ - w^—^ > — — Wp=n > 



- R=0.1 R=1.0 - R=3.0 — — R=0.1 

R=1.0 - R=3.0 


W C =\J .J wc— I -D Wc — J.J w c— V. J 

— We=l .5 - We=3.5 


- Nb=0.01 Nb=0.1 - Nb=0.3 - Nb=0.01 

— Nb=0.1 - Nb=0.3 


Figure 5: Thermal Radiation(R) 


Figure 6: Weissenberg Number (W) Figure 7: Effect of Brownian Motion 

(Nb) on Temperature in Different 
Electric Fields 


SWCNT s-water 

Solidline-El=0.0;Dashline-El=5.0 



- Nt=0.01 - Nt=0.1 - Nt=0.3 - Nt=0.01 

— Nt=0.3 - Nt=0.3 _ 


Figure 8: Thermophoresis 
Parameter (Nt) 


SWCNT s-water 

Solidline-El=0.0;Dashline-El=5.0 



n 


- Le=0.1 - Le=0.3 - Le=0.5 - Le=0.1 

■ ~ Le=0.3 - Le=0.5 _ 


Figure 9: Lewis Number (Le) 


SWCNTs-water 

Solidline-El=0.0;Dashline-El=5.0 



- X=0.1 - X=0.2 - a=0.3 - a=0.1 

■ — a=0.2 - a=0.3 


Figure 10: Effect of Chemical 
Reaction, on Temperature 
in Different Electric Fields 
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SWCNT s-\vater 

Solidline-El=0.0;Dashline-E 1=5.0 



- Y=0.1 - y=0.3 - y=0.5 - y=0.1 

- y=0.3 - y=0.5 _ 


Figure 11: Heat Generation (y) 
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SWCNT s-water 

Solidline-E 1 =0.0;Dashline-E 1=5.0 



- Pr=6.82 - Pr=9.82 - Pr=10..4 - Pr=6.82 

- Pt=9.82 - Pr=10-4 _ 


Figure 12: Effect of Parndtl Number 
(Pr) on Temperature in Different Electric Fields 


SWCNTs-water 

Solidline-E 1 =0.0;Dashline-E 1=5.0 



- a=0.1 - a=0.9 - X=1.5 - a=0.1 

— a=0.9 - a=1 .5 


Figure 13: Companion of the Chemical Reactions Obtained 
from Mair Khan et al. (2017) and the Current Study 

To confirm the present article with previous article, we compared the chemical reaction in articles which yalidates 
for this parameter. 


Table 2: Magnetic Field Deposition on / (0) 9 —0 (0) and —<p'( 0) 


M 

/"(0) 

-0'(O) 

-<P'( 0) 

Remark 

0 . 1 

-0.06192436 

0.974637508 

0.135958934 

E1 = 0.0 

0.5 

-0.21798656 

0.952379683 

0.151368698 


0.9 

-0.35423026 

0.933137979 

0.164942024 


0 . 1 

0.206321950 

1.014838665 

0.109158238 

E1 = 5.0 

0.5 

1.121600764 

1.139776354 

0.034226946 


0.9 

2.014618856 

1.244032672 

0.016290443 



Table 3: The Wall Thickness Deposition on / (0), —0 (0) and —<p\ 0) 


a 

/"(0) 

-0'(O) 

-<P'( 0) 

Remark 

0 . 1 

-0.06192436 

0.974637508 

0.135958934 

E1 = 0.0 

0.5 

0.027429696 

1.433762284 

0.151368698 


0.9 

-0.35423026 

0.933137979 

-0.25255788 


0 . 1 

0.098100963 

1.938636067 

-0.69151885 

E1 = 5.0 

0.5 

0.206321950 

1.014838665 

0.109158238 


0.9 

0.278513107 

1.468857871 

-0.27348653 
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Table 4: Power Low Index Deposition on / (0), —0 (0) and —<p'(0) 


n 

/"(0) 

-0'(O) 

-<p'( 0) 

Remark 

1. 2 

-0.06192436 

0.974637508 

0.135958934 

E1 = 0.0 

2 

-0.08075871 

0.970997376 

0.138390665 


3 

-0.35423026 

0.933137979 

-0.25255788 


1. 2 

-0.10151830 

0.967064068 

0.141038574 

E1 = 5.0 

2 

0.206321950 

1.014838665 

0.109158238 


3 

0.172513191 

1.008430946 

0.113222136 



Table 5: Thermal Radiation Deposition on / (0), —0 (0) and —<p'(0) 


R 

f"( 0) 

-0'(O) 

-<p'( 0) 

Remark 

0 . 1 

-0.06192436 

0.974637508 

0.135958934 

E1 = 0.0 

1.0 

-0.06192437 

0.833665065 

0.349542219 


3.0 

0.061924415 

0.716012830 

0.533825944 


0 . 1 

0.206321950 

1.014838665 

0.109158238 

E1 = 5.0 

1.0 

0.206321932 

0.863339187 

0.331765380 


3.0 

0.206321895 

0.733844123 

0.526822455 



Table 6: Weissenberg Number Deposition on /"(0), —0 (0) and —<p'(0) 


We 

f"( 0) 

-0'(O) 

-<P'(0) 

Remark 

0.5 

-0.08075871 

0.970997376 

0.138390665 

E1 = 0.0 

1.5 

-0.12009334 

0.963624344 

0.143371091 


3.5 

-0.18103546 

0.952941693 

0.150710643 


0.5 

0.172513191 

1.008430946 

0.113222136 

E1 = 5.0 

1.5 

0.111022612 

0.996790028 

0.120740391 


3.5 

0.026435463 

0.981229280 

0.131059540 



Table 7: Brownian Motion Deposition on / (0), —0 (0) and —<p'(0) 


Nb 

/"(0) 

-0'(O) 

-<P'( 0) 

Remark 

0.01 

-0.061924362 

1.027347931 

-5.33031215 

E1 = 0.0 

0.1 

-0.080758718 

0.970997376 

0.138390665 


0.3 

-0.061924371 

0.862442858 

0.540802231 


0.01 

0.2063219593 

1.068113220 

5.682418246 

E1 = 5.0 

0.1 

0.2063219507 

1.014838665 

0.109158238 


0.3 

0.2063219419 

0.901311372 

0.538079591 



Table 8: Thermophoresis Parameter Deposition on / (0), -0'(O) and -<p'(0) 


Nt 

f"(0) 

-0'(O) 

-<p'(0) 

Remark 

0.01 

-0.061924370 

1.017049917 

0.6766383438 

E1 = 0.0 

0.1 

-0.080758718 

0.970997376 

0.1383906658 


0.3 

-0.061924363 

0.890735511 

-1.036246600 


0.01 

0.2063219532 

1.112975344 

0.1006663215 

E1 = 5.0 

0.1 

0.2063219507 

1.014838665 

0.1091582385 


0.3 

0.2063219562 

0.931218226 

-1.143608929 



Table 9: Lewis Number Deposition on / (0), — 0'(O) and —<p'(0) 


Le 

/"(0) 

-0'(O) 

-<P'( 0) 

Remark 

0 . 1 

-0.061924368 

0.981158369 

0.0201155880 

E1 = 0.0 

0.3 

-0.080758719 

0.971571834 

0.1277775485 


0.5 

0.0619243645 

0.969535715 

0.2330978255 


0 . 1 

0.2063219488 

1.022196509 

-0.016019024 

E1 = 5.0 

0.3 

0.2063219505 

1.015493339 

0.0976142666 


0.5 

0.2063219538 

1.009105399 

0.2138938915 
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Table 10: Chemical Reaction Deposition onf (0), -0'(O) and -<p'(0) 


X 

f"( 0) 

-0'(O) 

-<p'( 0) 

Remark 

0 . 1 

-0.061924366 

0.974637508 

0.1359589349 

E1 = 0.0 

0.2 

-0.080758718 

0.969968474 

0.1698381885 


0.3 

-0.061924366 

0.972597445 

0.1982995950 


0 . 1 

0.2063219507 

1.014838665 

0.1091582385 

E1 = 5.0 

0.2 

0.2063219508 

1.013761455 

0.1407353665 


0.3 

0.2063219509 

1.012711055 

0.1717381518 



Table 11: Heat Generation Deposition on / (0), -0'(0) and -<p'(0) 


V 

f"( 0) 

-0'(O) 

-<p'( 0) 

Remark 

0 . 1 

-0.061924366 

0.974637508 

0.1359589349 

E1 = 0.0 

0.3 

-0.080758717 

0.710446099 

0.1653936047 


0.5 

-0.061924363 

0.396259884 

0.1956959508 


0 . 1 

0.2063219507 

1.014838665 

0.1091582385 

E1 = 5.0 

0.3 

0.2063219526 

0.768525605 

0.1304137086 


0.5 

0.2063219546 

0.470933108 

0.1562449437 



Table 12: Parndtl Number Deposition on / (0), -0'(0) and -<p'(0) 


Pr 

f"( 0) 

-0'(O) 

-<P'(0) 

Remark 

6. 82 

-0.06192436 

0.974637508 

0.135958934 

E1 = 0.0 

9.82 

-0.08075871 

1.153079098 

-0.13255348 


10.4 

-0.06192436 

1.191701399 

-0.18608932 


6. 82 

0.206321950 

1.014838665 

0.109158238 

E1 = 5.0 

9.82 

0.206321959 

1.206033800 

-0.16894764 


10.4 

0.206321960 

1.241286230 

-0.22020770 



4. CONCLUSIONS 

The present study was conducted to determine the effects of parameter such as thermal radiation, heat generation, 
thermophoresis, Brownian motion, Lewis number, magnetic field, wall thickness and power-law index on the heat transfer 
enhancement of nanoAuid water as base fluid. The numerical solution of the first order ODEs was obtained using the 
Runge-Kutta Fehlberg method with shooting technique. From the converged solutions of velocity, temperature and 
concentration the skin friction and heat and mass transfer rates were analyzed. The following conduction can be drawn. 


• The temperature decreases with respect to all parameters when the strength of the electric field is high (E 1=5.0) 
which means the electric field is dominate of water-SWCNT. 

• From the profiles of wall thickness, magnetic, radiation and heat generation strength, the thermal thickness is 
large. Therefore, its effect on temperature is quiet significant. 

• When electric field is strong, heat transfer increases with respect to parameters such as strength of magnetic field, 
Parandtl number and power-law index. Meanwhile heat transfer decreases with respect to parameters such as 
thermal radiation, Brownian motion, thermophoresis, heat generation, chemical reaction, Lewis and Wiessenberg 
number. 

• The results obtained from the preliminary tables show that the mass transfer rate increases with respect to electric 
field on thermal radiation, chemical reaction, heat generation, Lewis and Wiessenberg number. 
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